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Abstract. The Jeffrey-Kirwan residue is a powerful tool for computation of intersec- 
tion numbers or volume of symplectic quotients. In this article, we give an algorithm to 
compute it using Grobner bases. Our result is parallel to that of [CD] for Grothendieck 
residues. 

1. Introduction 

The Jeffrey-Kirwan residue was introduced in [JKilj and it is a powerful tool to 
compute intersection numbers or symplectic volume of symplectic quotients. There are 
several ways to compute it such as iterated residues, inverse Laplace transforms or nested 
sets f [JKi2] . [JKo] . [BV] . [SzV] . [dCP] 1. 

The aim of this article is to give an algorithm for computation Jeffrey-Kirwan residue 
using Grobner basis which can be implemented as a computer program. 

The contents of the article are as follows. In Section [2] we quickly review basic no- 
tions related to Grobner bases. In Section [3] we recall an algorithm for computation of 
Grothendieck residue from |CD] which is similar to ours. In Section H] we generalize ideas 
from previous section which will be the core of our algorithm. In Section [5] we recall the 
definition and properties of Jeffrey-Kirwan residue. In Section [6] we apply the results of 
Section [2] in the case of Jeffrey-Kirwan residue and we give an algorithm to compute it 
using Grobner bases. 

Acknowledgement. The author is grateful to Andras Szenes and Michele Vergne 
for discussions and useful comments. The support of FNS grant 132873 is gratefully 
acknowledged. 

2. Grobner bases and the division algorithm 

In this section we recall some notions and results from the book [GLOj about Grobner 
bases and division of multivariate polynomials. Let IK be a field (R or C) and consider 
the polynomial ring K[x\, ... ,x r ]. We use the notation x a , a = (a±, . . . , a r ) G N r for the 
monomial x® 1 ■ ■ ■ , as special case x° = 1. 

Definition 1. A monomial order " > " on . . . ,x r ] is an order on monomials with 

properties: 

(i) it is a total order, 

(ii) x a > x b implies that x a+c > x b+c for all c £ N r , 

l 
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(iii) x a > 1 for all a £ N r \{0}. 

Example. The lexicographic order is defined as follows. Let a = (01, . . . ,a r ) and b = 
(pi, ... , b r ) in N r . Then x a >i cx x b if in a — b = (a± — b±, . . . , a r — b r ) the left-most 
non-zero entry is positive. 

Definition 2. Fix a monomial order and let / = X^aeA c a xa be a non-zero polynomial 
in K[a;i, . . . , x r \. Let x m be the greatest monomial of the set {x a \ a £ A} with respect to 
the fixed monomial order. Lm(/) = x' m is called the leading monomial of /. Moreover, 
we call Lt(/) = c m x m the leading term of / and Lc(/) = c m is called the leading 
coefficient of /. Remark that Lt(/) = Lc(/) • Lm(/). If / is an ideal in K[xi, . . . ,x r ] 
then the ideal of leading terms Lt(7) = (Lt(/) | / € I) is the ideal generated by the 
leading terms of all polynomials in /. 

Definition 3. A finite subset G of / is called a Grobner bases of I with respect to the 
chosen monomial order if Lt(/) = (Lt(<?) | g £ G). 

It is not unique since any finite subset of / containing a Grobner bases is also a Grobner 
bases. However, for each polynomial ideal with a monomial order there is a unique 
reduced Grobner bases G of I with properties: 

(i) Lc(#) = 1 for all g £ G, 

(ii) for all g £ G, no monomial of g belongs to (Lt(G \ {g}))- 

Grobner bases of an ideal in . . . , x r ] can be computed by Buchberger's algorithm, 

see [UL"Uj. §11.7. 

Next, we recall the division of multivariate polynomials, see |GLOj . §11.3. We fix 
a monomial order onK[ii,...,x r ] and let F = {fi, ■ ■ ■ , f s } be an ordered s-tuple of 
polynomials. Using the following division algorithm any polynomial / can be written as 
/ = ai/i + . . . + a s f s + R with a\, . . . , a s , R € K[x\, . . . , x r ] such that no monomials of 
R is divisible by any of Lt(/i), . . . , Lt(/ s ). The polynomial R is called the remainder of 
the division of / by F. The division algorithm goes as follows. In the beginning we set 
R to f and a±, . . . , a s to 0. Suppose we have written 

/ = ai/i + • • • + a s f s + R 

for some oi, . . . , a s , R. Choose the greatest monomial M in R (with respect to the fixed 
monomial order) such that it is divisible by Lt(/j) with smallest i. Then set the new ctj 
to dj + j^hry and the new R to R— jjfijjfi- The algorithm ends in finite steps since the 
order of monomials in R which is divisible by any of Lt(/i), . . . , Lt(/ s ) strictly decrease 
in each step and monomial orders are also well-orders. 

The remainder of the division in general depends on the way the elements of F are 
listed. Moreover, it is possible that the remainder is non-zero even if / £ (/i, . . . , / s ). 
An example is given in Example 5 of |CLOj . §11.3. Nevertheless, if F = {/i, . . . , f s } is 
a Grobner bases of the ideal • • • , / s ) then the remainder R does not depend on the 
way the elements of F are listed. 
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Definition 4. Fix a monomial order. The normal form Nj(f) of / with respect to the 
ideal / is a polynomial such that / — Nj(f) € / and no monomial of Nj(f) is contained 
in Lt(J). 

Proposition 5. IfG = {g±, . . . ,g s } is a Grobner bases of the ideal I, then the remainder 
of the division of f by G is equal to the normal form of f , i.e. R = Np[f). Moreover, 
f £ I if and only if R = (compare |CLUj . §11.6, Proposition 1 and Corollary 2). 

Remark. If / is a homogeneous polynomial ideal, then for any monomial order we can 
choose a Grobner bases G = {g±, . . . ,g s } with gi homogeneous polynomials. Moreover, 
the remainder of the division of a homogeneous polynomial / by G has is also homo- 
geneous of the same degree as /. Furthermore, Ni(f) is also homogeneous of the same 
degree as /. 

Lemma 6. Let I be a homogeneous polynomial ideal with a monomial order. Then I 
and Lt(J) are isomorphic as graded vectors spaces, i.e. dim/^ = dimLr(/)d for all d 
(more generally, |CLOj . §IX.3, Proposition 4). 

3. Review of Grothendieck residue case 

In this section we review an algorithm for Grothendieck residue as explained in [CD], 
§1.5.6. Let H, P\, . . . , P r E C[xi, . . . , x r ] be polynomials of degree d, d±, . . . , d r , respec- 
tively. 

Definition 7. Suppose that G C r is the only common zero of P±, . . . , P r in an open 
neighborhood U of it. The Grothendieck residue is defined as 
/ H \ 1 f H{x) 

ReS ° \P,\...\Pr) ~ (&7=Tf Me) Pl(x) ■ ■ ■ P r (x) ?1 • • • 

where r(e) = {x £ U | |Pj(a;)| = £i > 0, i = 1, . . . ,r} and it is oriented by (argPi(x)) A 
. . . A ( arg P r (x)). 

The Grothendieck residue has the following properties. 

(i) It is independent on e. 

(ii) If P%, . . . , P r and H are homogeneous, using rescaling x i— > Xx we can see that 

Res o(prrW) = 0ifd ^EI=i(^-i)- 

(iii) (Local duality) If Reso ( p^rrp ) = for all G € C[x±, . . . ,x r ] then H lies in the 
ideal {P 1 ,...,P r ). 

(iv) (Transformation law) Let Qi, ■ ■ ■ ,Q r £ C[xi, . . . ,x r ] such that is their only 
common zero locally. If we can write Qi = Y^j=i a ijPj then 

(v) (Euler- Jacobi vanishing theorem) Let Pi, . . . ,P r such that G C r is the only 
common zero. If d < Y^i=i(di ~ 1) then Reso ( \p ) = 0- 
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As noted in |CD] the Grothendieck residue can be computed using normal forms as 
follows. Suppose that Pi, . . . , P r are homogeneous with only common zero £ C. 

(1) Compute a Grobner bases of the ideal I = (Pi, . . . , P r ). 

(2) Compute normal forms Nj(H) of H and Nj(A) of A = det j=i) , where 
JPj = y^-_-| dijXj, (i = 1, . . . , r). 

(3) Reso( Pi | g |p ) = jVj(A) ' shice by transformation law Reso( Pl | A |p ) = Reso( 3 ,i 1 < x ) 
= 1. 

The non-homogeneous case can be reduced to the homogeneous case as follows. Let 
Pi, . . . ,p r G C[xi, . . . , x r ] be polynomials of degree di, . . . , d r such that their only com- 
mon zero is G C r . Let 

Pi(x ,xi, . . . ,x r ) = x$ ■ Pi (—,...,— ] G C[x ,a;i, . . . ,x r ) 

\Xq XqJ 

be the homogenization of pi. Let h G C[zi, . . . ,x r ] be a degree d polynomial. If d < 
YJi=i(di — 1) then Reso(^j^-) = by Euler-Jacobi theorem. Therefore, suppose that 
d > YH=i(di — 1) and let 

H(x ,xi, ... ,x r ) = Xq ■ h (—,..., — 

be the homogenization of h. Moreover, let Pq(xq,x\ . . . ,x T ) = Xq° with do = d + 1 
ELi(^-l)- Then 

h \ ( H 

.pJ^Wr) ~ ReS ° Ul^l-I^ 

and we can apply the above algorithm to compute it. 

4. The main idea 

Let I be a field and denote the polynomial ring . . . ,x r ]. Denote IfCfx]^ 

the set of homogeneous polynomials of degree d. Moreover, let A : — > K be a 
non-zero linear functional. We extend it to the whole polynomial ring by zero. 
Let I = {P G K[x] | X(P ■ Q) = 0, \/ Q G K[x]} which is a homogeneous polynomial 
ideal. Remark that it is enough to check if X(P ■ Q) = only for monomials Q. To 
show that / is homogeneous, suppose that P £ I and decompose it to homogeneous 
components P n + P n -\ + . . . + Pq, and for all Q = Q m homogeneous of degree m we 
have A(P • Q) = \{Pd-m ■ Qm) = 0. For degree reasons, A(P^-m ■ Qm 1 ) — for all Q m i 
homogenenous of degree m' ^ m. Hence Pd-m G I and all homogeneous components of 
P are contained in /. In particular, 1^ the degree d part of I is equal to the kernel of A 
in K[x] d . 

Proposition 8. Let A G K[x] d such that A(A) / 0. Then Nj(A) / and for all 

P G K[x]d we have 

A(P) NijP) 
A(A) JVj(A)" 
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Proof. Since A(A) 7^ 0, A is not contained in / and by Proposition [5j Nj(A) is a non- 
zero homogeneous polynomial of degree d. Since Id is one codimensional in ICfxj^, the 
vector space LT(/)d is also one codimensional. Therefore, for all P e X^J^ the normal 
form Nj(P) is a multiple of the same monomial M of degree d which is not contained in 
LT(/) d , i.e. N^P) = Dr(iV/(P)) = Lc(JV/(P)) • M. Finally, for all P e K[x] d we have 

X(P) = \{Nj(P)) = Lc(JV/(P)) • A(LM(iV/(P))) = hc(Ni(P)) • A(M), 

and 

A(P) _ Lc(AT 7 (P)) • A(M) _ Lc(iV/(P)) • M _ jVj(P) 
A (A) ~ Lc(iV/(A)) • A(M) ~ Lc(iV/(A)) • M ~ Nj(A)' 

□ 

If we can compute the ideal / and a non-zero A(A) for some A then we can also compute 
A(P) for all P G Xfa;]^ using the following algorithm: 

(1) compute a Grobner bases G of /, 

(2) compute normal forms Nj(P) and Nj(A) as remainders of the division of P and 
A by G, 

(3) finally, A(P) = A(A) • G X. 

In the case of Jeffrey-Kirwan residue we have to compute the ideal I and we have to also 
find a A. 

5. The Jeffrey-Kirwan residue 

We recall the definition and properties of Jeffrey-Kirwan residue following | JKoj . Let 
V be an r-dimensional real vector space and let 21 = [a±, . . . , a n ] be a collection of (not 
necessarily distinct) non-zero vectors in V*. We consider a^s as linear functions on V. 
Let A be a connected component of V\U2 =1 a^~, where af- = {v G V \ a%{v) = 0}. Remark 
that for all i, either G A v or -Qj G A v , where A v = {/3 G V* \ (3(v) > 0, We A} 
is the dual cone of A. If aj e A v then we say that it is polarized. Choose ( £ A and 
a bases {x%, . . . ,x r } of V* such that = 1, and X2(£) = ■■■ = x r (£) = 0. Let 

e = £\Xi + . . . + e r x r e V* and let P e R[V] be a polynomial. We define 

' P(x)e e ^ 
if ei < 0, 



+ I P(x)e^ 

Milium)- 1 



Res Xl=00 — 

— -x x if £1 > 0, 



JKRes ( " v '~ x ) = = Res+ ( . . . ( Res+ 

\n i= i<xi(x) j Jdet[( Xi , Xj )}? j=1 v v Ui=i»i(x) t 



considering X2, ■ ■ ■ ,x r as constants while taking the residue with respect to x\. Fix a 
scalar product on V* and define the Jeffrey-Kirwan residue 

P(x)e £ ^ \ _ 1 + ( (^ + P{x)e £ ^ 

Definition 9. We call an element of V* regular if it does not lie on any (r — 1)- or less 
dimensional subspace of V* spanned by a subset of 21. 
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The Jeffrey-Kirwan residue has the following properties. 

(1) Suppose that a±, . . . , a r G A v and suppose that e is regular. Let P G R[V] be a 
homogeneous polynomial of degree d. Then 

JKRes A ( -ffi^ x) = 

unless all the following properties are satisfied 

(a) {«i, . . . , a n } span V* as vector space, 

(b) d < n — r, 

(c) e G Cone(ai, . . . , a n ). 

(2) Suppose that d <n — r and P is homogeneous of degree d. Then 



JKRes A 



P(x)e- 



e(x) 



xWlim JKRes^flW^; 

7 ^^°+ v fc! nr=i«i(^) , 



JKRes A 



e(x) n - r - d P(x)e £ ^ 



(n-r-d)\ Ui=i a i( x Y 
(3) If d = 0, re = r and properties (l)(a) — (c) are satisfied then 

e*(aO \ i 



JKRes 



A 



ni=iai(^)7 I det(«i, . . . ,a r 



where [a±, . . . , a r ) denotes the matrix whose columns are the coordinates of 
ai, . . . ,a r with respect to any orthonormal bases of V*. 

Remark. For e regular the Jeffrey-Kirwan residue does not depend on the choice of £ in 
A and the corresponding bases {x\, . . . , x r } of V*. 

We turn to a more algebraic approach used in |BVj . Suppose that we have fixed a 
polarization A and a±, . . . ,a n are polarized with respect to it. Let i?a = R|X]si be the 
set of rational functions of form p - , where P G WL\V] and m,- > 0. 

Definition 10. A subset / C {1, . . . ,n} is called generating if {ai \ i £ 1} span V as 
vector space. Moreover, a subset I C {1, . . . , n} is called e-generating if it is generating 
and e £ Coneicii \ i G I). 

We have the following partial fraction decomposition (more generally, Theorem 1 of 

Mb- 

Proposition 11. Let G% C Pg be spanned as vector space by fractions of form 1 mi 

lli=l a i 

such that the set I = {i \ rrii > 0} is generating. Moreover, let NG% C R% be the 

vector space spanned by fractions of form p m- such that the set I = {i \ mi > 0} is 

1 l 1= i a i 

non- generating. We have a direct sum decomposition 

Pa = G a © NG*. 
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Denote (R%)d the degree d part of R%, d G Z. In particular, (-Ra)d — (G'a)d © (-W^sOd 
and (Gg)_ r is spanned by fractions 1 such that {ai-,, . . . , ctj } is a bases of V*. In 

general, these fractions are not linearly independent, for example 



1 



x(x+y) xy (x+y)y~ 

Definition 12. Fix a scalar product on V*. Let e £ V* be regular and consider the 
linear functional JK e on R% which vanishes on R<& \ {G%)- r and on (Ga)-r is defined by 

1 \ I 7-777 " 77 if e G Cone(a h ,...,a ir ), 

JK £ ( I = < |det(a il) ...,Q!i P )| 

otherwise, 



OLj r i ' ' ' OLj, 



and (aiu • • • , cti r ) denotes the matrix whose columns are the coordinates of . . . , ai r 
with respect to any orthonormal bases of V*. Moreover, we introduce a shorter notation 
-21/ d\ _ m ( P 



JKf(P) = JK £ ^ 
Remark. 

(1) If P € K[V] is a homogeneous polynomial of degree n — r then 

p(a;) e e(aO \ P 



JK Res A AT j ; , x = JK 



and for other polynomials is related via property (2) of the Jeffrey-Kirwan 
residue. 

(2) In the case of V* = W with the standard scalar product, 21 = [ai, . . . , a n ] a gen- 
erating subset of Z r and e regular the Jeffrey-Kirwan residues can be interpreted 
as intersection numbers on the toric variety X c (2l): 

JK e (y^—) = [ X(P), 

\lU=l a iJ JX C (%) 

where c is the connected component of regular elements containing e and \ '■ 
] — > iJ*(X c (2t)) is a degree preserving ring homomorphism (degXj = 
2) as explained in |SzV| . §1 and §2. 

6. Computation of Jeffrey-Kirwan residue using Grobner bases 

To apply the result of Section 2] in the case of Jeffrey-Kirwan residue, we have to 
compute the ideal I %£ = {P G R[V] \ JK*(P ■ Q) = 0, VQ G R[V]} and we have to 
compute a non-zero JK^(A) for some A G R|V]. Since e G Cone(a.i \ i = 1, . . . , n), 
there are a^, . . . ,aii r such that e G Cone(a{ 1 , . . . , aii r ). By the property (3) of the 
Jeffrey-Kirwan residue JK e f- — ^— — ) = 1-5—7 — tt 7^ 0, therefore let A = . 

Theorem 13. We have = (lljeJ a i I £ ^ Cone(aj \ j ^ J)). More geometrically, 
the ideal is generated by all ria 4 eH+ a * suc ^ that H is a hyperplane spanned by some 
of ai 's and H + denotes the open half-space containing e. 
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Proof. Introduce the notation I^ £ = {Y[j^j a j \ £ & Cone(aj \ j ^ J)). The inclusion 
1% £ C igi,e is straightforward: if P = ]X g j ctj is a generator of e then for all Q we 
have JK £ ( n £^ a . ) = JK^^-^-) = 0, since e Ctme(aj | j £ J). 

We show the inclusion igt,e C e by induction on n. If n = r then there are two 
cases: 

(1) if {1, . . . , n} is not e-generating then = WL[V] = 1% £ , 

(2) if {1, . . . , n} is e-generating, then igt )£ = (a±, ... ,a r ) = I% £ . 

Suppose that n > r and {1, . . . ,n} is e-generating. Denote 2lj := 21 \ and remark 
that there is an index i such that 2lj is e-generating. We can think of igi i)£ = i ^ 



R[F] | JK £ (jg^) = JK £ (n^u) = 0, VQ G R[V]}, hence I a , £ C J a . £ for all 
i = 1, . . . , n. By the induction hypothesis 

where 1%. £ = ( Yij^j oij \ i ^ J and {1, . . . , n} \ ({i} U J) is not e-generating). We will 
show that 

Define the homomorphism it : M.[yi, . . . ,y n ] — > M[V], 7r(yj) = Oj. Introduce the ideals 
^St,e = ( IljeJ yj \ £ ^ Cone(otj \ j ^ J)) and similarly define Lgi ii£ for all i = 1, . . . , n. 
The ideals L$ £ and e are square- free monomial ideals (i.e. the generators are square- 
free monomials) and n^ =1 Lgi i£ is also square-free monomial ideal. Moreover, Tr(L% t£ ) = 
Iq £ and 7t(L<h u£ ) = J^. £ . Therefore, it is enough to show that 

L%, e = n™ =1 L*2 lij£ . 

Take a square-free generator Hfee-ft" °^ ^^=i-^sti,e> ^ C {1, . . . , ra}. Remark that the 
index set K ^ § since there is z such that 2lj is e-generating. Fix an m E K. Since 
Ylk£K Vk £ L% mtE , there is a subset J m C {1, . . . , n} such that m ^ J m and TljeJ m G 
L*2 lmj£ divides Ilfee-ftr^' Hence J m U {m} C -ftT, but the complement of J m U {m} 
in {l,...,n} is not e-generating, implying that the complement of K is also not e- 
generating, i.e. U keK yk G This shows n™ =1 La ij£ C L« )£ . 

To show that Lg[ )£ C r\f =1 L% ii£ , let HjeJ^i e ^2t,e- Then the complement of J is not 
e-generating. Let 

k J t £ J ' 

then rijej, % e ^a i)£ which divides nj 6 j%> nence ELfe,/^ G for a11 i - 

To prove the second part of the theorem, denote H% t e the ideal generated by Jl Qi g/f+ on 
such that H is a hyperplane spanned by some of a^s and H + denotes the open half- 
space containing e. It is easy to see that -ff*a,£ C igt,e- To show the reverse inclusion let 
rijejOj G 7gj )£ , i.e. J C {1, . . . , n} is a subset such that e ^ Cone(aj \ j ^ J). We can 
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suppose that J is minimal, hence the latter convex polyhedral cone is at least (r — 1)- 
dimensional. Therefore, it has an (r — l)-dimensional face such that the underlying 
hyperplane H separates e from this convex polyhedral cone. Moreover, H is spanned 
by some aj's and {ctj \ j ^ J} C {ctj \ ctj ^ H + }, i.e. {ctj \ ctj G H + } C {ay \ j G J}, 
hence ria!jeff+ a i divides njeJ a i> therefore n^gj^i e ^2t,e- D 

Remark. The ideal 7<a,e is the image of the Stanley-Reisner ideal Iy, of a fan £ under 
the map 7r : R[yi, . . . , y n ] — > ^[x±, . . . , x r ], j/i i — >• Oj. To construct the fan £ consider the 
short exact sequence 

^ M. n - r W 1 — W 

where 7 sends the standard basis element ej to Qj. Its Gale dual 

^ W R n — ^ R n ~ r ^ 

comes with vectors /3\, . . . ,/3 n G K n_r as /3j = <5*(ej). The fan £ is given as follows: for 
a subset J C {1, . . . ,n\ the Cone(f3j \ j G J) G £ if and only if e G Coneicti \ i ^ J). 
Since e is regular £ is a simplicial fan. To the fan £ we associate the Stanley-Reisner 
ideal / s = < Vj I Cone(/3j \ j G J) £ £> C . . . , y„]. 

Example. In the r = 2 case Me is a hyperplane of y* and suppose that vectors a±, . . . , 
and «fc+i; • • • > a n he in different components of l/*\Me. Then I% e = (a\ ■ ■ ■ a^, ctk+i ■ ■ ■ ct n ) 

Finally, the Proposition [5] and Theorem [T3] with the remark preceding it gives the 
following proposition about computation of Jeffrey-Kirwan residues. 

Proposition 14. Let 21 = [ct\, . . . , a n ] be a collection of non-zero vectors in V* , dim V* = 

r. Let A be a connected component of V \ Uf =1 a^, suppose that ct\, . . . ,a n G A v and 
e G V* is regular. Let J C {1, . . . , n} be a subset such that {ctj \ j G J} is a bases of V* 
with e G Coneipij \ j G J). T/ien 



ni=i a i0)7 ' |det( aj |j G J)| iV/^dl^j «i) 

/or a// P G M[V] homogeneous of degree n — r and (ctj \ j G J) denotes the matrix whose 
columns are the coordinates of ctj , j G J with respect to any orthonormal bases of V* . 

Moreover, JKRes I y-,„ ' — i-^-x can be computed with the following algorithm: 
\il;=i a i\ x ) ' / 

(1) Compute the ideal L% j£ using Theorem [T51 

(2) Find a subset J C {1,... ,n} such that {ctj \ j G J} is a bases of V* with 
e G Cone{a.j \ j G J). 

(3) Compute a Grobner bases G of igt,e- 

(4) Compute the remainders (normal forms) Ni sle (P) and Ni m e (Y\ i( £j ct-i) of the 
division of P and Yl^ j cti by G, respectively. 
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(5) JKRes 



,A 



( 



nr=i a i( x ) ' 




) 



det(oy \j e J)\ N I%e (TJ^j on) ' 



1 N hl jP) 



Remark. The remainders can be computed by computer programs like Maple or Macaulay2 
as follows. Let / be an ideal of R[a;i, . . . , x r \. In Maple, first compute the Grobner bases 



with(Groebner) : G=Basis(I,tdeg(x_l, . . . ,x_r)) 
Then compute the remainder of the division of / by G: 
NormalForm(f ,G,tdeg(x_l , . . . ,x_r) ) 

In Macaulay2 we can compute the remainder directly, 
f % I 
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